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In  this  paper,  we  develop  tractable  mathematical  models  and  approximate  solution  algorithms 
for  a  class  of  integer  optimization  problems  with  probabilistic  and  deterministic  constraints,  with 
applications  to  the  design  of  distributed  sensor  networks  that  have  limited  connectivity.  For  a 
given  deployment  region  size,  we  calculate  the  Pareto  frontier  of  the  sensor  network  utility  at  the 
desired  probabilities  for  d-connectivity  and  fc-coverage.  As  a  result  of  our  analysis,  we  determine 
(i)  the  number  of  sensors  of  different  types  to  deploy  from  a  sensor  pool,  which  offers  a  cost 
vs.  performance  trade-off  for  each  type  of  sensor,  (ii)  the  minimum  required  radio  transmission 
ranges  of  the  sensors  to  ensure  connectivity,  and  (iii)  the  lifetime  of  the  sensor  network.  For  gen¬ 
erality,  we  consider  randomly  deployed  sensor  networks  and  formulate  constrained  optimization 
techniques  to  obtain  the  localization  performance.  The  approach  is  guided  and  validated  using 
an  unattended  acoustic  sensor  network  design.  Finally,  approximations  of  the  complete  statis¬ 
tical  characterization  of  the  acoustic  sensor  networks  are  given,  which  enable  average  network 
performance  predictions  of  any  combination  of  acoustic  sensors. 

Categories  and  Subject  Descriptors:  C.2.1  [Computer-Communication  Networks] :  Distributed 
networks;  G.1.6  [Numerical  Analysis]:  Optimization — Constrained  optimization,  convex  pro¬ 
gramming,  integer  programming,  nonlinear  programming ;  G.3  [Probability  and  Statistics]: 
Experimental  design 

General  Terms:  Algorithm,  Design,  Performance. 

Additional  Key  Words  and  Phrases:  Bayesian  experimental  design,  dynamic  programming,  sensor 
networks. 


1.  INTRODUCTION 

The  design  and  operation  of  sensor  networks  are  rich  with  integer  and  combinatorial 
optimization  problems.  Some  of  these  problems  concern  (i)  the  dynamic  assignment 
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of  a  subset  of  already  deployed  sensors  for  target  tracking  to  minimize  sensor  net¬ 
work  power  consumption,  [Liu  et  al.  20031  ISubhlok  et  al.  19991  ISinha  and  Chandrakasan  2001], 
(ii)  the  automatic  placement  of  a  given  set  of  sensors  to  guarantee  coverage  and 
connectivity,  [Cowan  and  Kovesi  19881  llsler  et  al.  2004j.  and  (iii)  the  optimal  move¬ 
ment  strategies  for  mobile  sensors  on  a  graph  in  the  context  of  probabilistic  pursuit 
evasion  games,  [Vidal  et  al.  2002:  Isler  et  al.  2004|.  Despite  the  advances  that  have 
been  made,  many  difficult  issues  remain  to  be  solved.  The  common  theme  of  these 
sensor  network  problems  is  that  the  available  resources  are  already  spent  on  a  set 
of  sensors.  Consequently,  it  is  desired  to  obtain  the  best  knowledge  about  a  state- 
of-nature  with  minimum  effort  given  the  available  sensors. 

In  this  paper,  we  provide  statistical  models  and  mathematical  programming  so¬ 
lutions  to  determine  the  optimal  sensor  choices  for  deployment  from  a  given  sensor 
pool  under  a  limited  budget.  Therefore,  we  focus  the  available  resources  on  a  set  of 
complimentary  sensors  for  the  deployment  purpose  before  the  sensors  are  deployed. 

Optimality  is  obtained  by  simultaneously  maximizing  three  criteria  under  a  desired 
minimum  connectivity  level  of  the  sensors.  These  criteria  are  (1)  a  desired  utility 
of  the  sensor  network  deployment,  (2)  a  desired  minimum  lifetime  of  the  sensor 
network,  and  (3)  a  desired  sensor  coverage  on  the  deployment  area.  For  generality, 
we  assume  random  sensor  deployment,  where  each  sensor’s  position  comes  from 
a  homogenous  Poisson  point  process  over  an  area.  We  refer  to  the  multi-criteria 
optimization  as  the  network  design  strategy  (NDS)  problem. 

In  |  Cevher  and  Kaplan  2008] ,  we  develop  a  theory  to  predict  the  localization 
performance  of  randomly  distributed  sensor  networks  consisting  of  various  sensor 
modalities  when  only  a  constant  active  subset  of  sensors  that  minimize  localization 
error  is  used  for  estimation.  The  characteristics  of  the  modalities  include  mea¬ 
surement  type  (bearing  or  range)  and  error,  sensor  reliability,  FOV,  sensing  range, 
and  mobility.  We  show  that  the  localization  performance  of  a  sensor  network  is  a 
function  of  a  weighted  sum  of  the  total  number  of  each  sensor  modality.  We  also 
show  that  the  optimization  of  this  weighted  sum  is  independent  of  how  the  sensor 
management  strategy  chooses  the  active  sensors. 


In  this  paper,  we  define  a  sensor  network’s  localization  utility  in  a  Bayesian  ex¬ 
perimental  design  framework  [Chaloner  and  Verdinelli  1995).  We  then  employ  our 
results  in  determining  the  localization  performance  of  heterogenous  sensor  networks 
in  Cevher  and  Kaplan  2008']  to  justify  the  form  of  the  localization  utility.  We  then 
discuss  and  summarize  how  the  sensor  network’s  lifetime  depends  on  the  connec¬ 
tivity  of  the  sensors  and  what  coverage  implies  for  collective  parameter  estimation. 
Different  from  j  Cevher  and  Kaplan  2008] ,  we  also  provide  continuous  relaxation  so¬ 
lutions  that  alleviate  the  computation  and  provide  insight  to  the  final  form  of  the 
NDS.  While  the  focus  of  (Cevher  and  Kaplan  2008]  is  the  derivation  of  the  localiza¬ 
tion  accuracy  of  sensor  networks,  this  paper  focuses  on  the  design  aspects  of  sensor 
networks  and  how  they  can  be  efficiently  obtained. 


Our  approach  is  guided  and  validated  by  the  design  of  acoustic  sensor  networks 
connected  with  wireless  links.  Acoustic  sensor  networks  use  single  microphone 
sensors  and  tethered  acoustic  arrays  to  detect,  track,  and  classify  acoustically  loud 
targets  and  events  such  as  ground  vehicles,  helicopters,  and  sniper  fire.  Target 
detection  and  parameter  estimation  performance  of  the  acoustic  sensors  increase  as 
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the  number  of  microphones  in  the  sensor  is  increased  (Johnson  and  Dudgeon  1 993  i 
ICevher  2005] .  Hence,  it  is  possible  to  create  a  sensor  pool  using  acoustic  sensors 
with  different  number  of  microphones. 

Our  specific  contributions  are  as  follows: 


(1) 

(2) 


We  provide  approximations  to  the  statistical  characterization  of  acoustic  sen¬ 
sors.  These  approximations  enable  average  performance  predictions  for  sen¬ 
sor  networks  that  consist  of  any  combination  of  acoustic  sensors.  In  contrast, 

we  use  generic  sensing  models. 


Cevher  and  Kaplan  2008 


We  provide  analytical  solutions  based  on  branch-and-bound  relaxations  that 
can  upper  and  lower  bound  the  integer  programming  solutions.  These  approx¬ 
imations  enable  quick  calculations  of  optimal  sensor  network  designs. 


(3)  We  provide  theorems  that  underlie  the  elements  of  the  optimal  sensor  network 
design:  concentration  of  resources  and  dominating  sensor  pairs.  These  condi¬ 
tions  reduce  the  search  space  for  the  optimization  problems  and  enable  fast 
calculations.  Specifically,  we  prove  that  given  a  total  of  Nu  objectives,  the  op¬ 
timal  sensor  network  design  chooses  at  most  Nu  different  sensor  types  in  the 
solution.  We  also  prove  that  under  conditions  described  in  Section  [3731  there 
exists  sensor  pairs  which  automatically  eliminate  other  sensors  in  the  Pareto 
frontier. 


(4)  We  extend  the  dynamical  programming  formulation  in  Cevher  and  Kaplan  2008| 
to  account  for  multiple  utilities. 


The  organization  of  the  paper  is  as  follows.  In  Sect.O  we  describe  the  NDS  prob¬ 
lem  in  the  context  of  Bayesian  experimental  design,  d-connectivity,  and  k- coverage 
problems.  In  Sect.  [3J  we  describe  optimization  solutions  to  the  NDS  problem.  In 
Sect.lU  we  derive  the  performance  metrics  for  acoustic  sensor  networks.  In  Sect.0 
we  demonstrate  the  NDS  solution  for  acoustic  sensor  networks  for  position  estima¬ 
tion  where  we  determine  the  Pareto  frontier  for  utility  and  lifetime  of  the  sensor 
network. 


2.  THE  SENSOR  NETWORK  DESIGN  OBJECTIVES 

The  sensor  network  design  simultaneously  optimizes  multiple  objectives  that  sum¬ 
marize  various  performance  criteria  related  to  wireless  sensor  networks.  These 
objectives  are  formulated  as  a  function  of  the  characteristics  and  quantities  of  each 
feasible  sensor  type  and  how  the  sensor  network  is  managed.  Typical  characteris¬ 
tics  of  the  sensor  types  include  measurement  modality  (bearing  or  range)  and  error, 
sensor  reliability,  field-of-view,  sensing  range,  and  mobility.  The  design  output  is  a 
vector  n  =  [ni  ri2  •  •  ■  tit}'  ,  whose  elements  consist  of  the  number  of  each  sensor 
type  from  a  catalogue  of  T-sensors,  denoted  as  s  =  {st\t  =  1, . . .  ,T},  for  deploy¬ 
ment.  Each  sensor  type  has  a  different  monetary  cost  and  only  a  finite  budget  is 
available  to  purchase  and  operate  these  sensors.  The  sensor  costs  are  expressed  as  a 
T-dimensional  vector  c,  whose  t-th  element  Ct  corresponds  to  the  cost  of  St  .  Hence, 
a  valid  sensor  design  must  satisfy  the  budget  constraint  c'n  <  $. 

The  sensor  network  design  relies  on  the  ability  to  use  the  characteristics  of  the 
sensor  types  to  predict  the  localization  accuracy,  lifetime,  coverage,  and  reliability 
of  the  sensor  network  for  any  given  combination  of  sensor  types  and  deployments 
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satisfying  the  budget  constraint.  This  mapping  involves  three  entangled  aspects 
of  sensor  networks:  deployment  plans,  operational  choices,  and  sensor  choices.  In 
order  to  make  the  problem  tractable,  we  assume  that  the  sensors  are  deployed  ran¬ 
domly  over  a  coverage  area  A  of  size  A  via  a  uniform  distribution.  By  ignoring 
boundary  effects,  we  assume  that  the  node  positions  are  realizations  of  a  homoge¬ 
nous  Poisson  point  process  (PPP)  [Ross  2007] .  For  the  operation  of  the  sensor 
network,  a  sensor  management  strategy  based  on  selecting  a  total  of  g-sensors  to 
minimize  the  localization  error  is  assumed.  Furthermore,  the  t-th  sensor  type  has 
an  effective  sensing  range  of  Rf,  which  is  important  in  obtaining  expressions  for 
k- coverage.  Finally,  the  transmission  of  data  is  the  dominant  user  of  battery  re¬ 
sources,  and  the  transmission  range  rtran  is  determined  based  upon  the  density  of 
the  nodes  to  ensure  a  desired  d-connectivity  of  the  network.  It  should  be  noted 
that  the  objectives  (or  utilities)  are  defined  to  be  non-negative  real  numbers  where 
larger  numbers  represent  higher  performance. 

2.1  Localization  Accuracy 

The  purpose  of  the  a  sensor  network  is  to  exploit  the  disparate  data  collected  from 
the  various  nodes  in  order  to  better  infer  the  state  of  the  objects  (or  targets)  9 
comprising  a  scene,  e.g.,  position,  velocity,  appearance,  physical  structure,  pitch 
frequency,  loudness,  etc.  Before  the  actual  deployment  of  the  sensor  network,  tar¬ 
gets  are  assumed  to  be  or  to  appear  within  the  area  A,  and  their  states  have  a  prior 
distribution  p(9),  e.g.,  targets  have  certain  speed  distributions,  time-frequency  and 
loudness  characteristics.  After  deployment,  the  sensor  network  provides  noisy  ob¬ 
servations  of  9  via  the  sensory  outputs,  making  it  possible  to  judge  the  amount  of 
information  provided  by  the  sensor  network  over  the  prior  knowledge  of  9.  In  this 
paper,  our  primary  focus  is  the  localization  accuracy  of  the  sensor  network  on  A. 

After  deployment,  the  location  of  the  sensor  nodes  can  be  represented  as  the  set 
C  =  {Ct|i  =  1,  ■■■,?},  where  Ct  =  [  Ct,i  Ct,2  ■  ■  •  C t,nt  ]  is  the  matrix  of  the  loca¬ 
tions  of  all  nt  sensors  of  type-f  and  C  is  the  2 -D  vector  of  actual  sensor  locations. 
Each  sensor  provides  data  y  =  {yk,t\k  =  1, ...  ,71*;  t  =  1, . . . ,  T},  e.g.,  acoustic  mi¬ 
crophone  recordings  within  some  dynamic  range.  The  characteristics  of  the  sensors 
determine  a  probability  density  function  (PDF)  of  the  measurements  as  a  function 
of  the  sensor/target  geometry,  i.e.,  p(y\C,9.n).  Then,  the  measurements  can  be 
used  to  provide  an  estimate  of  the  target  state  9.  Note  that  9  can  be  viewed  a 
function  of  the  measurements  y,  where  the  shape  of  the  function  is  dictated  by  the 
measurement  PDF  p(y|C,  9 ,  n).  As  the  data  is  not  observed  before  n  is  chosen  and 
an  arbitrarily  large  number  of  real-world  realizations  is  possible,  it  is  meaningful 
to  choose  n  to  minimize  the  average  estimation  error  of  the  target  state  over  all 
random  realizations 

e(n)  =  J  J  J  p{e(y)\C,6,n)p(6)p(t)e(yA,O,n)dyd0d£,  (1) 

where  e(y,  9:  n)  represents  the  accuracy  of  the  state  inference: 

e(y,C,  9,n)  =  (e-d(y^j  A  (e -0(y)^  ,  (2) 

and  A  is  a  positive  definite  weighting  matrix.  In  this  work,  we  consider  the  mean 
squared  error  (mse)  so  that  A  =  J,  where  I  is  the  identity  matrix. 
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When  the  relationship  between  y,  9,  n,  and  £  is  nonlinear,  approximations 
based  on  Taylor  series  expansions  are  typically  used  to  avoid  the  complicated  inte¬ 
gral  in  dT|)  |Berger  1993}  IChaloner  and  Verdinelli  19951  Tierney  and  Kadane  1986 


These  approximations  use  the  expected  Fisher  information  matrix  or  the  matrix  of 
the  second  derivatives  of  the  log  likelihood  function  and  are  extremely  accurate  even 
for  small  data  sample  sizes.  Hence,  assuming  (i)  a  network  design  that  guarantees 
v  >  3  sensors  contributing  to  the  estimation  at  any  spatial  location  (fc-coverage), 
(ii)  a  network  with  high  connectivity  (d-connectivity),  and  (iii)  that  each  sensor  has 
sufficiently  large  number  of  data  samples,  we  can  use  the  following  normal  approx¬ 
imation  for  the  data  likelihood  in  <JT])  (AT  (mean,  variance)  is  the  Gaussian  density): 


p(0{y)\0X,n)=N(e,  [H  +  3{0,C,n)\  ')  ,  (3) 

where  9  is  the  mode  of  the  posterior  distribution,  H  is  the  Hessian  of  the  loga¬ 
rithm  of  the  prior  density  p(9)  (precision  matrix  of  the  prior),  and  $(9,  £,  n)  is  the 
expected  Fisher  information  matrix  (FIM)  IChaloner  and  Verdinelli  1995j.  When 
the  prior  for  the  target  state  is  uninformative,  e.g.,  uniform,  H  =  0.  While  some 
domain  knowledge  such  as  a  road  may  lead  to  the  use  of  an  informative  prior, 
this  paper  assumes  that  no  such  knowledge  is  available:  H  =  0  and  p{6)  =  1/A. 
With  the  normal  approximation  l(3|).  the  expected  utility  0  can  be  expressed  as 
follows  (Chaloner  and  Verdinelli  1995] : 

e(")  JJ  tr{5_1(0,C,n)}p(C)d0dC-  (4) 

We  define  the  localization  utility  as  the  reciprocal  of  the  expected  MSE  as  defined 
in  0.  Typically,  the  MSE  expression  is  too  complex  to  compute  in  closed  form. 
Because  the  targets  and  sensor  nodes  are  distributed  uniformly  over  A.  we  model 
the  sensor/target  geometry  using  2 -D  PPP.  As  a  result,  the  actual  MSE  scales  with 
the  density  of  the  network  normalized  to  the  quality  of  the  sensor  measurements. 
As  shown  in  |Cevher  and  Kaplan  2008;  Kaplan~and  Cevher  2007|,  when  the  sensors 
are  stationary,  the  localization  utility  (within  a  linear  scale  factor)  is  given  by: 

ftnt]  ,  (5) 


U(n)  =  (£ 
\t= 1 


where  ft  depends  on  the  inherent  reference  sensing  accuracy  at,  field-of-view  at  € 
(0, 1],  and  reliability  probability  /3t  €  (0, 1]  of  the  sensor  type-t: 


ft 


at  Pt 
2/7 ' 


(6) 


The  derivations  of  0  and  ©  assume  that  only  g-active  sensors  are  used  for  esti¬ 
mation  at  each  snapshot.  The  scale  factor  of  0  depends  on  the  number  of  active 
sensors  q  and  what  objective  the  sensor  management  optimizes  in  active  sensor 
selection.  In  Cevher  and  Kaplan  2008],  utility  expressions  including  mobile  sen¬ 
sors  are  also  provided.  The  parameter  7  >  0  depends  on  the  signal  propagation 
loss  factor  of  the  medium.  For  the  case  of  no  sensor  management  where  all  nodes 
within  the  sensing  range  to  the  target  participate  in  the  localization  given  that  the 
obtain  at  detections,  the  MSE  is  approximately  inversely  related  to  0.  Monte 
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Carlo  experiments  in  Section  [5]  demonstrate  that  the  reciprocal  of  ©  is  indeed  a 
good  approximation  to  the  MSE.  Therefore,  we  use  m  as  the  localization  accuracy 
utility.  The  specific  weights  for  acoustic  sensors  are  obtained  computationally  in 
Section  0 

2.2  Effects  of  Connectivity  and  Lifetime 

Connectivity  (C)  is  an  important  property  of  multi-hop  wireless  networks,  where 
each  sensor  cooperates  in  routing  each  other’s  packets.  In  distributed  sensor  net¬ 
work  operation,  it  is  desirable  that  the  sensors  do  not  communicate  raw  data  and 
the  communication  among  the  sensors  has  a  constant  bandwidth  to  maximize  the 
sensor  lifetime  and  to  minimize  interference  among  the  receivers  sharing  the  same 
wireless  channel  when  the  sensors  are  communicating  or  when  the  network  is  queried 
for  information  |  Akyildiz  et  at.  2002| .  In  addition,  the  network  estimates  are  based 
on  local  sensor  observations,  which  can  determine  the  parameter  of  interest  6  after 
fusion.  For  example,  when  bearing  and  range  sensors  are  used,  determining  the 
target  position  relies  on  collective  observability,  which  requires  connectivity. 

In  graph  theory  terms,  the  network  is  d-connected  (with  connectivity  degree 
C  =  d)  if,  for  any  given  pair  of  sensor  nodes,  there  exists  at  least  d  mutually 
independent  paths  connecting  them  (Bollobas  19981.  Each  sensor  usually  has  a 
fixed  transmission  range  rtran >  which  is  significantly  smaller  than  the  dimensions  of 
A.  A  communications  link  can  be  created  between  two  sensors  only  if  their  physical 
distance  is  less  than  their  transmission  range  (Fig.QJa)).  Hence,  given  transmission 
range  rtran)  the  probability  of  creating  links  between  sensors  increase  as  the  number 
of  sensors  in  the  field  increases. 


Fig.  1.  (a)  A  sensor  node  (vertices)  can  connect  to  its  neighbors  by  creating  links  (edges)  only  if 
other  nodes  are  within  its  transmission  range  rtran-  Then,  to  connect  to  any  other  node  in  the 
network,  the  sensor  can  multi-hop  its  information  using  its  neighboring  nodes,  (b)  Approximate 
connectivity  probability  vs.  the  number  of  sensors  for  a  sensor  network  q  with  the  following 
parameters:  A  =  lOOOmx  1000m  and  rtran  =  100m. 


Under  the  random  deployment  assumption,  it  is  possible  to  show  that  the  pdf 
of  the  nearest  neighbor  sensor  distance  p  is  given  by  p(p)  =  2nXpexp  (— 7rAp2) , 
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where  A  =  ^  is  the  spatial  node  density,  and  N  =  ]Ct=i  nt  (total  number  of 
sensors)  [Bettstetter  2002].  Hence,  the  probability  that  each  sensor  has  at  least 
d-neighbors  within  the  transmission  range  rtran  (node  degree  D  =  d)  on  a  toroidal 
surface  is  given  by  the  following  [Bettstetter  2002] : 

P{D>  d\n,()  =  exp  <|  N  log  ^1  -  ^ exp  (-7rArt1 2ran)  j  |  .  (7) 

Figure  [0(b)  gives  example  plots  of  the  node  degree.  It  is  shown  in  [Penrose  1999| 
that  in  a  graph,  as  the  number  of  vertices  increase,  the  node  degree  D  converges  to 
the  connectivity  degree  C  with  probability  one.  In  this  paper,  we  approximate  the 
connectivity  of  the  sensor  network  by  0.  In  general,  the  connectivity  probability 
should  be  simulated  for  the  specific  geometry  of  A  (e.g.,  see  [Bettstetter  2002)  for 
simulations  on  square  *4).  It  is  important  to  note  that  the  connectivity  probability 
monotonically  increases  as  a  function  of  N.  Hence,  to  maximize  the  connectivity 
probability  at  any  given  degree,  the  total  number  of  sensors  N  in  the  sensor  network 
must  be  maximized. 

Impacts  of  Connectivity  on  Power :  Note  that  the  connectivity  can  always  be 
increased  by  increasing  the  sensor  transmission  range  rtran.  However,  the  choice  of 
v tran  also  affects  the  lifetime  (L)  of  the  sensor  network^]  A  sensor  consumes  power 
for  (,4.)  transmitting  data,  (B)  receiving  data,  (C)  sensing,  (X>)  aggregating/fusing 
information,  (£)  idling,  and  [IF)  coping  with  radio  interference/communications 


overhead/etc.  [Heinzelman  et  al.  2002|  Bhardwaj  et  al.  2001: 

Zhang  and  Hou  2004|| 

Bhardwaj  and  Chandrakasan  2002;  Feeney  and  Nilsson  2001 

Lindsey  et  al.  2002  . 

The  sensor  transmission  range  rtran  directly  affects  A.  and  also  indirectly  affects  IF. 


In  |Bhardwaj  et  al.  2001| ,  A-C  are  assumed  to  derive  geometry-independent  tight- 
bounds  for  the  lifetime  of  a  sensor  network,  where  L  oc  IV.  The  reason  for  this 
proportionality  is  that  the  total  energy  of  the  sensor  network  scales  with  the  total 
number  of  sensors.  If  the  communication  protocols  are  smart  enough  to  minimize 
losses,  then  the  actual  lifetime  should  follow  the  same  trend.  The  bound  is  not 
tight  when  the  number  of  sensors  is  small  as  the  characteristic  distance,  needed 
to  make  the  cost  of  transmitting  a  bit  linear  with  distance,  cannot  be  realized. 

In  [Bhardwaj  and  Chandrakasan  2002| ,  D  is  also  assumed,  and  an  improved  bound 
is  given  which  is  also  geometry  dependent.  In  [Zhang  and  Hou  2004|,  A-£  is  as¬ 
sumed  for  randomly  deployed  sensor  networks  and  it  is  shown  that  the  lifetime  is 
also  almost  linear  with  the  total  number  of  sensors  when  algorithms  that  can  adap¬ 
tively  turn  on/off  sensors  for  energy  conservation  are  used.  In  [Giridhar  and  Kumar  2005| , 
it  is  shown  that  message  passing  among  the  nearest  neighbors  yields  nearly  optimal 
lifetime. 

The  explicit  dependence  of  the  sensor  network  lifetime  on  rtran  and  N  is  beyond 
the  scope  of  this  paper.  This  is  because  the  analysis  requires  an  underlying  network 
model,  which  can  quantify  the  MAC  collision  and  communication  delay  issues. 
Without  focusing  on  the  specifics  of  the  network  models,  we  argue  that  the  typical 
MAC  collision  and  communication  delay  issues  have  negligible  effect  on  the  sensor 
network  lifetime  in  our  case  because  we  only  use  q  <C  N  active  sensors  per  target 


1  The  lifetime  can  be  defined  in  various  ways  such  as  the  time  to  first  node  failure  due  to  battery 

depletion  or  the  time  to  appearance  of  the  first  connectivity  brake-down. 
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for  localization  at  any  given  time.  Moreover,  we  assume  that  the  sensor  network 
queried  at  random  sensor  locations,  that  is,  there  is  no  fixed  base  station;  hence, 
the  traffic  load  and  the  energy  consumption  per  sensor  are  uniformly  distributed 
across  the  sensor  network.  Then,  based  on  the  results  from  the  literature,  we  can 
expect  have  a  monotonic  behavior  for  L  as  a  function  of  N . 

In  the  acoustic  sensor  network  design,  when  pp  =  1  -  e  (e  <  1),  we  assume  that 
the  maximum  network  lifetime  has  approximately  the  following  form: 

LocN6,  (8) 

where  6  >  0  for  generality  and  the  proportionality  is  independent  of  rtran.  This 
form  is  motivated  by  the  idealized  case  where  the  sensor  batteries  are  depleted  only 
by  the  /'^-propagation  loss  ( p  >  4):  if  pc  is  high,  then  L  oc  C~1(N)  x  r~p ,  where 
C(N)  is  the  average  hop  count.  For  a  square  A,  the  average  hop  count  increases 
approximately  proportional  to  N _1/2  with  N,  whereas  r~p  oc  Np/2  (see  footnote  2). 
Hence,  L  oc  N^p~1^2 .  Other  forms  for  L  can  be  assumed,  which  result  in  tractable 
solutions,  e.g.,  L  oc  exp  (SN).  In  general,  we  emphasize  again  that  the  lifetime  has 
a  monotonic  (quasiconvex)  form  for  any  given  pc  for  efficient  solutions. 


Fig.  2.  (a)  Connectivity  probability  pc  as  a  function  of  N  and  rtran-  (b)  Sensor  network  lifetime 
L  as  a  function  of  N  and  rtran. 

Figures  [2]( a)  and  (b)  roughly  illustrate  the  interplay  between  the  sensor  network 
lifetime  and  connectivity.  As  the  total  number  of  sensors  is  increased,  the  desired 
connectivity  probability  pc  can  be  achieved  with  smaller  sensor  transmission  ranges 
Ftran-  For  a  given  transmission  range  that  allows  high  connectivity,  the  lifetime 
of  the  sensor  network  linearly  increases  with  the  total  number  of  sensors  in  the 
network  as  the  total  energy  of  the  network  is  linearly  increased  (assuming  efficient 
communication  protocols  that  can  minimize  losses).  Since  the  connectivity  can  be 
preserved  with  smaller  transmission  ranges  as  the  number  of  sensors  is  increased@ 


2  As  an  example,  if  pc  =  1  —  e  (e  <g  1),  then  from  0  with  d  —  1,  WG  have  Ttran  — 

\j - — — ^ ^ — -  Si  —  | A\  -log  oc  TV-1/2  for  | log  e |  3>  | log  JV | ,  which  is  satisfied  for 

typical  N. 
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the  actual  lifetime  of  the  sensor  network  grows  faster  than  the  growth  obtained  by 
increasing  N  while  keeping  rtran  fixed. 

2.3  Effects  of  Coverage 

The  coverage  ( V )  problem  in  sensor  networks  aims  to  quantify  how  well  the  area 
of  interest  A  is  monitored.  The  coverage  problem  has  been  extensively  studied  in 
the  literature,  see  Sect.  2  in  [Lazos  and  Poovendran  2006]  for  a  survey.  Given  a 
planar  A  and  random  sensor  network  deployment,  one  can  determine  the  number 
of  sensors  so  that  any  given  point  in  A  is  sensed  by  at  least  fc-sensors  with  a  given 
probability  (^-coverage).  For  the  sensor  network  design,  the  k- coverage  probability 
is  quite  complicated  as  each  sensor  can  have  a  heterogeneous  sensing  range. 

In  [Altman  and  Miorandi  2005j.  a  Poisson  approximation  is  given  for  the  k- 
coverage  probability  for  randomly  deployed  sensor  networks  with  heterogeneous 
sensing  ranges: 

P(V  >  k\n ,  C)  =  — ,  where  v  =  ^  Ktnt ,  Kt  =  (9) 

l>k  '  t=  1 

Note  that  the  coverage  probability  is  a  monotonically  increasing  function  of  v.  To 
increase  k- coverage  probability  to  a  desired  level  py,  v  in  ([9])  must  be  maximized. 

3.  OPTIMIZATION  PROBLEMS 

The  goal  of  NDS  is  to  maximize  three  competing  objectives:  1)  localization  accu¬ 
racy,  2)  lifetime,  and  3)  coverage.  The  exact  form  of  these  three  objectives  were 
discussed  in  the  previous  section.  NDS  is  basically  a  multi  criteria  optimization 
problem  over  the  simplex  A f  defining  the  feasible  set  of  network  designs,  i.e.,  the 
vectors  n  whose  elements  are  non-negative  integers  and  also  satisfy  the  budget  con¬ 
straint  c'n  <  $.  Because  of  the  integer  nature  of  the  network  design  vectors  n,  the 
number  of  elements  in  the  simplex  A f  is  finite,  albeit  a  very  large  number.  In  gen¬ 
eral,  the  objectives  (or  utilities)  are  expressed  as  an  Nu  dimensional  vector  u  whose 
elements  represent  the  Nu  different  utility  values.  For  NDS,  Nu  =  3.  Furthermore, 
rti,  zt2,  and  U3  represent  the  localization  accuracy,  lifetime,  and  coverage  utilities, 
respectively. 

The  following  subsections  describe  general  multiple-objective  optimization,  exact 
integer  programming  methods  to  find  the  “optimal”  designs,  and  approximate  non¬ 
integer  optimization  methods  to  quickly  determine  a  “nearly  optimal”  design. 

3.1  Multi-Objective  Optimization 

Multi-objective  problems  are  typically  addressed  using  Pareto  optimality.  In  other 
words,  an  optimal  network  design  is  one  such  that  one  cannot  find  another  sensor 
design  in  A f  where  one  or  more  of  the  objectives  increase  but  none  decrease.  The 
locus  of  all  possible  optimal  network  designs  forms  the  Pareto  frontier,  i.e,  the  trade¬ 
off  surface  for  the  utilities.  In  other  words,  the  Pareto  frontier  is  the  collection  of 
feasible  utilities  that  cannot  be  uniformly  dominated  by  any  other  feasible  utility. 
In  general,  enumeration  of  all  possible  parameter  choices  that  satisfy  the  constraints 
is  required  to  explore  the  Pareto  frontier.  However,  in  the  end,  some  preference 
must  be  made  to  choose  an  operating  point. 
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Many  points  on  the  Pareto  frontier  can  be  determined  via  scalarization  (or  weight¬ 
ing),  where  one  simply  maximizes  a  weighted  sum  of  the  utilities  ATu  for  any  A 
whose  elements  are  non-negative  |Boyd  and  Vandenberghe  2004j.  The  scalarization 
concepts  also  hold  for  the  maximization  of  a  product  of  utilities 


Nu 

IK.  (10) 

i= 1 

because  the  log  operator  is  monotonically  increasing.  Figure  [3]  illustrates  an  exam¬ 
ple  Pareto  surface  for  two  objectives.  If  the  Pareto  frontier  happens  to  be  concave, 
then  all  the  points  can  be  determined  via  scalarization. 


Fig.  3.  An  example  efficient  frontier  between  log  U  and  logL  is  shown.  Shaded  region  is  the 
feasible  region  for  the  problem  that  can  be  achieved  by  increasing  rtran-  The  dashed  part  of  the 
Pareto  frontier  cannot  be  determined  using  scalarization  as  there  is  no  supporting  plane  to  the 
set  of  feasible  U  and  L  values.  Among  the  supporting  planes,  A'  maximizes  U  X  L. 


We  consider  four  general  optimization  problems  relevant  to  the  NDS: 

-Problem  1:  Determine  the  Pareto  frontier  for  the  Nu  objectives, 

-Problem  2:  Maximize  one  objective  over  M  without  regard  to  the  values  of  the 
other  Nu  —  1  objectives, 

-Problem  3:  Maximize  the  product  scalarization  of  the  utilities  for  a  given  A, 
and 

-Problem  4:  Maximize  one  objective  while  constraining  the  other  objectives  to 
exceed  a  threshold. 

Once  Problem  1  is  solved,  it  is  straightforward  to  solve  the  other  three  single 
objective  problems.  For  Problems  2-3,  one  simply  determines  the  point  on  the 
Pareto  frontier  that  maximizes  the  scalar  quantity  of  interest.  Finally,  Problem  4 
searches  for  the  maximal  point  on  a  subset  of  the  Pareto  frontier  where  the  con¬ 
straints  are  satisfied. 

In  general  the  relationship  between  the  utilities  and  the  network  design  is  a 
vector  of  nonlinear  functions,  i.e. ,  u  =  h(ni, . . .  ,um)-  As  shown  in  Section [2J  all 
three  NDS  objectives  can  be  expressed  as  a  monotonically  increasing  function  of  a 
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weighted  sum  of  sensor  type  populations  for  a  given  network  design,  i.e. , 


Ui  = 


for  i  =  1, . . . ,  3. 


Note  that  the  weights  u>itt  are  all  non-negative. 


(11) 


3.2  Integer-Programming 

In  this  section,  we  discuss  a  dynamic  programming  solution  to  enumerate  the  pa¬ 
rameter  values  comprising  the  Pareto  frontier,  thereby  solving  Problem  1.  Note 
that  the  Pareto  frontier  is  a  subset  of  the  surface  of  the  feasible  region  (see  Figure[3| 
where  the  feasible  region  is  the  locus  of  utility  values  corresponding  to  all  sensor 
designs  in  the  simplex  A/”.  It  is  useful  to  tabulate  the  associated  network  design  and 
cost  to  a  given  utility  in  the  feasible  region.  Therefore,  we  define  the  feasible  triple 
as  the  set  of  triplets  describing  the  utility,  NDS  vector,  and  cost  for  each  point  in 
the  feasible  region  as  follows 

T  =  {(utility:  h(n),  design:  n,cost:  c'n )  :  n  €  A/”}.  (12) 

Clearly,  the  feasible  triple  as  given  by  (H2D  can  be  generated  by  considering  each 
network  design  in  A/",  and  the  feasible  region  is  simply  Tu  =  ^.utility,  i.e., 

Au  =  {{.utility  :  f  £  T}.  (13) 

The  finite  cardinality  of  the  simplex  M  leads  to  a  finite  cardinality  for  !FU  and  T . 

Fortunately,  the  relationship  between  the  network  designs  and  the  utilities  is 
given  by  m.  and  one  need  not  resort  to  an  exhaustive  search  over  the  simplex 
A f.  Many  of  the  triples  in  the  feasible  region  are  redundant  in  the  sense  that 
they  can  be  removed  in  T  without  affecting  the  feasible  region  J~u  given  by  m 
Formally,  a  triple  fi  £  is  redundant  if  there  exist  another  triple  {2  such  that 
f  1  .utility  =  {2  .utility  but  fi.cost  >  f2.cost.  The  observation  that  the  design  ni 
leads  to  a  redundant  triplet  (h(ni),  ni,  c'ni),  i.e.,  there  exist  another  design  71.2 
such  that  h(rii)  =  h(ri2)  but  dn-i  <  c'ni,  instantly  distinguishes  many  other 
network  designs  as  also  redundant.  For  instance,  the  network  design  n  1  +  A,  where 
the  elements  of  A  are  non-negative  integers,  leads  to  a  triplet  that  is  redundant  in 
light  of  the  triplet  associated  to  the  design  n2  + A.  Algorithm[T]is  a  dynamic  integer 
program  to  determine  the  feasible  region  that  exploits  the  redundancy  property  to 
prune  the  search  space. 

In  Algorithm  [1]  0  is  a  column  vector  of  zeros  and  et  is  the  unit  vector  whose 
f-th  element  is  one  and  whose  other  elements  are  zero.  Parameter  Nf,max  is  the 
feasible  design  with  the  most  number  of  sensors,  and  n  represents  the  total  num¬ 
ber  of  sensors  associated  to  the  network  designs  being  analyzed.  The  sets  Co  and 
Ci  represents  the  set  of  triples  associated  to  network  designs  of  n  —  1  and  n  sen¬ 
sors,  respectively,  that  are  currently  non-redundant,  and  the  set  C  accumulates  the 
non-redundant  triples.  For  each  non-redundant  design  consisting  of  n  sensors,  the 
algorithm  sequentially  adds  one  sensor  of  type  t  to  the  design  in  line  6  and  tests  if 
such  a  design  leads  to  a  redundant  utility  vector  via  lines  7  and  8.  If  the  design  is 
currently  non-redundant  and  meets  the  cost  constraint,  its  triple  is  added  to  Ci  and 
C.  Furthermore,  the  redundant  triple  is  removed  from  C.  Note  that  the  redundancy 
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test  is  performed  before  applying  the  nonlinearity  in  CD-  Therefore,  C  accumu¬ 
lates  the  minimal  set  of  network  designs  to  generate  the  feasible  set,  and  line  17 
transforms  this  set  into  the  minimal  feasible  triple  by  incorporating  the  nonlinear 
function  h(-). 


Algorithm  1. 


(1)  Determine  -/V/-jlnax 


J  where  c*  =  mint  Ct  ■ 


(2)  Initialize  C0  =  {(Ojv.xi,  0Txi,  0)}, C  =  0,Ci  =  0. 

(3)  for  n  =  1, . . . ,  max, 

(4)  for  t  =  l,...,T, 

(5)  for  all  f*  £  Co, 

(6)  f+  =  (f*. utility  +  flet,  f*.  design  +  et,  f.cost  +  ct) ; 

(7)  if  3f  £  C  such  that  f+. utility  =  f. utility, 

(8)  if  f+. cost  <  f.cost , 

(9)  remove  f  from  C  and  insert  f+  into  C  and  C\ ; 

(10)  elseif  f+.cost<$, 

(11)  insert  novel  f+  into  C  and  C\ ; 

(12)  end; 

(13)  end; 

(14)  Co  <—  Ci ; 

(15)  Ci=0; 

(16)  end; 

(17)  T  =  {(h(f.utility),  (.design,  f.cost)  :  f  £  C}  forms  the  minimal  triple. 


Algorithm  [T]  requires  that  the  relationship  between  the  utilities  and  the  network 
design  is  given  by  m ,  but  it  does  not  assume  any  more  structure  regarding  the 
utilities.  For  the  NDS,  the  second  element  of  the  utility  vector  is  the  lifetime 
objective,  which  is  equivalent  to  the  total  number  of  sensors.  Thus,  the  second 
element  of  the  utility  vectors  associated  to  Co  and  C\  are  always  n  —  1  and  n 
respectively.  The  test  in  line  7  can  only  be  positive  if  f  £  C\. 

Before  the  application  of  Algorithm  [I]  the  weighting  elements  LUi}t  should  be 
transformed  into  integers  by  an  appropriate  scaling  within  a  desired  level  of  accu¬ 
racy.  Integer  w^t’s  in  turn  lead  to  integer  values  for  (.utility.  Furthermore,  these 
values  can  be  bounded  by  iV&  =  Nf,max  max<x>ijt.  In  other  words,  the  localization 
accuracy  and  coverage  utilities  form  unique  values  and  the  existence  and  non¬ 
existence  of  these  Njf  unique  utility  values  in  Cq  and  C\  can  be  stored  in  AT&  x  AT& 
matrices.  This  replaces  the  implicit  search  in  line  7  with  a  table  look-up.  As  a 
result,  the  overall  complexity  of  Algorithm [l]  is  0(Nf^maxTN^). 

Once  the  non-redundant  feasible  triple  is  determined  by  Algorithm[lJ  it  is  straight¬ 
forward  to  determine  the  Pareto  frontier.  The  Pareto  frontier  V  is  simply  the  max¬ 
imal  subset  of  such  for  any  u  £  V. utility  and  any  vd  £  futility  with  u/u', 
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then  it  is  not  the  case  that  zq  <  u\  for  all  i  6  {1 , ,M},  i.e. ,  u  is  not  pointwise 
less  than  u ' .  A  simple  method  to  compute  V  is  to  exhaustively  remove  any  point 
in  iFu  that  is  pointwise  less  than  another.  Let  Vi  represent  the  surface  formed  by 
maximizing  the  i-th  objective  over  points  in  J~n  whose  other  objectives  are  equal. 
In  other  words,  if  u  £  7^. utility,  €  iFu,  and  Uj  =  uj  for  j  ^  i,  then  Ui  >  u\. 
Clearly,  if  one  arbitrarily  selects  a  utility  u  from  T. utility  there  exist  an  utility  ir 
from  Vi- utility  that  is  pointwise  greater  than  or  equal  to  u.  Thus,  Vt  is  a  superset  of 
V.  As  demonstrated  in  Section  0  this  feasible  surface  is  not  necessarily  equivalent 
to  the  Pareto  frontier.  As  mentioned  at  the  end  of  Section  EU  once  the  collection 
of  points  forming  the  Pareto  frontier  is  known  (Problem  1),  then  the  solutions  to 
Problems  2-4  are  straightforward. 

3.3  Non-integer  Optimization 

The  optimal  NDS  solution  to  Problem  1  determines  the  Pareto  frontier  where 
feasible  designs  occupy  a  simplex  A/”  in  the  T-dimensional  integer  lattice  ZT.  By 
relaxing  the  simplex  J\f  in  the  T-dimensional  real  space  RT,  optimization  algorithms 
can  become  more  computationally  efficient.  The  solutions  obtained  by  continuous 
relaxations  are  usually  not  optimal.  However,  when  the  total  budget  $  is  rather 
large,  optimal  solutions  select  the  key  sensor  types  in  great  numbers.  Then,  the 
difference  between  the  solutions  obtained  by  continuous  relaxations  and  the  optimal 
integer  solutions  become  negligible.  Continuous  relaxations  also  allow  us  to  gain 
further  insights  into  the  other  NDS  problems,  as  we  will  see  this  section. 

For  the  ease  of  notation,  we  define  the  following  transformations:  nm  =  cmnm/$ , 
&i,t  —  Uij/Cm,  and  [f2]jt  =  wq*  so  that  the  NDS  problem  is  equivalent  to 

Find  the  Pareto  frontier  of  f In  subject  to  l'h  <  1  and  h  >  0,  (14) 

where  1  is  a  T  dimensional  column  vector  of  ones.  The  following  three  theorems 
alleviate  the  search  for  the  optimal  points  on  the  Pareto  frontier. 

Theorem  1 .  ( Concentration  of  Resources )  The  optimal  NDS  design  n*  asso¬ 
ciated  with  the  Pareto  frontier  of  the  multi-objective  optimization  problem  in  0 
expense  the  available  resources  in  at  most  Nu  different  sensor  types. 

Proof.  The  Pareto  frontier  is  a  subset  of  feasible  utilities  in  the  positive  orthant 
of  the  7Vu-dimensional  space.  Within  the  feasible  region,  any  point  can  be  viewed 
as  a  unit  vector  times  a  scale  factor.  A  point  in  the  feasible  region  where  the  scale 
factor  is  maximized  for  a  given  direction  is  located  on  the  Pareto  frontier  of  (fT4l). 
Let  us  consider  the  an  arbitrary  unit  vector  in  the  positive  orthant,  i.e.,  e  such  that 
e  y  0  and  e'e  =  1.  Then,  consider  the  following  the  linear  program 

min  l'n  subject  to  fin  =  e  and  h  >;  0.  (15) 

This  linear  program  has  Nu  equality  constraints;  hence,  the  optimal  solution  h* 
contains  at  most  Nu  nonzero  elements  (Theorem  20.1)  [Moon  and  Sterling  2000], 
Let  $*  =  l'n*  be  the  minimum  objective  value  of  the  linear  program  in  (1151)  at  the 
optimal  parameter  value  n*.  Then,  a  network  design  sh*  where  s  =  l/$*  meets  the 
cost  constraint  of  m  and  leads  to  a  feasible  point  se  that  lies  on  the  boundary. 
Otherwise,  if  there  was  another  design  fi\  that  leads  to  a  feasible  point  in  the  same 
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direction  but  larger  scale  factor,  then  fi*  could  not  solve  the  linear  program  in  (USD- 
Namely,  if  Ori-i  =  sie  where  Si  >  s  and  Vfii  <  1,  the  NDS  design  n 2  =  l/sifti 
would  solve  the  constraint  in  (usd  at  a  cost  of 

l'fi2  =  l/sil'ni 

<  1/si, 

<  1/s  =  $*  =  1'h* , 

which  creates  a  contradiction  because  h*  is  the  optimal  solution  of  m- 

The  proof  is  true  for  any  direction  in  the  positive  ortlrant;  hence,  any  network 
design  associated  with  the  Pareto  frontier  has  at  most  Nu  nonzero  elements.  This 
result  reduces  the  NDS  sensor  search  space  from  2T  =  Y^t=  1  (I)  to  Nsc  =  (  ^  ) 
sensor  combinations.  □ 

Theorem  2.  (Dominating  Sensor  Pairs-1)  Consider  two  indices  ( j,k )  £  [1,T] 
and  j  <  k  for  which  the  row  vectors  forming  the  matrix  Cl  in  G2P  are  locally 
strictly  convexWm  £  ( j,k ),  i.e.,  WijTTl  <  Qij  +  Then,  no  sensor  m 

is  assigned  any  resources  on  the  Pareto  frontier.  That  is,  sensors  j  and  k  dominate 
any  sensor  m  in  the  solution. 

Proof.  Assume  a  feasible  n*  for  which  h ^  >  0  and  define  a  new  resource 
distribution  where  hj  =  h*  +  (dh^,  hk  =  n%  +  (1  —  /3)n^,  and  hm  =  0  while  keeping 
all  other  elements  the  same  from  n* ,  where  (3  =  .  The  utilities  associated  to 

these  two  designs  are  u*  =  Tin*  and  u  =  Sin,  and  (m  —  u*)  ln*m  =  POjij  +  (1  — 
0)u>itk  ^  &i,m-  Thus  by  the  convexity,  u,  >  u*  for  all  objectives  i.  Hence,  on  the 
Pareto  frontier  =  0,  otherwise  it  is  always  possible  to  increase  all  the  objectives 
in  the  multi-criteria  optimization  problem. 

If  the  dominating  sensor  pairs  are  on  the  boundary  of  the  sensor  book,  then  the 
final  solution  of  the  sensor  network  design  has  at  most  three  configuration  choices: 
only  sensor  type  t  =  1,  only  sensor  type  t  =  T,  or  a  combination  of  sensor  types 
t  =  1  and  t  =  T.  □ 

Theorem  3.  (Dominating  Sensor  Pairs-2)  Consider  two  indices  ( j,k )  £  [1,T] 
and  j  <  k  for  which  Nu  —  1  row  vectors  forming  the  matrix  Tl  in  are  locally 
strictly  convex  Vm  £  ( j ,  k)  with  respect  to  the  ordering  of  remaining  row  vector, 
i.e.,  uji^m  <  z’k  u>i  j  +  tbiki^i  7^  i' ■  Then,  no  sensor  m  is  assigned 

any  resources  on  the  Pareto  frontier.  That  is,  sensors  j  and  k  dominate  any  sensor 
m  in  the  solution. 

Proof.  The  proof  is  similar  to  that  of  Theorem[2j  and  it  is  omitted.  □ 

The  difference  between  Theorem  2  and  Theorem  3  is  the  ordering  of  the  feasible 
sensors.  Whereas  the  ordering  of  the  sensors  in  Theorem  2  is  arbitrary,  Theorem 
3  requires  that  they  are  sorted  with  respect  to  the  order  of  one  of  the  rows  of 
Tl.  Theorems  1-3  provide  the  following  simplifications  of  the  NDS  problems  in 
Section  l3Tl 

3.3.1  Solution  of  Problem  2.  In  light  of  TheoremQ]  the  solution  of  Problem  2 
simply  assigns  all  the  resources  to  a  single  sensor  type  because  Nu  =  l,  i.e.,  h  =  et 
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for  some  t.  The  sensor  modality  that  maximizes  the  utility  is  given  by 

t*  =  argmaxuq^.  (16) 

Hence,  one  simply  chooses  the  modality  that  provides  the  biggest  bang  for  the  buck. 


3.3.2  Solution  of  Problem  3.  The  solution  to  Problem  3  also  assigns  resources 
to  at  most  Nu  sensor  types  due  to  Theorem  1.  Theorems  2  and  3  can  be  used  to 
further  decrease  Nsc  by  removing  sensors  that  are  dominated  in  the  feasible  sensor 
set.  Subsequently,  given  each  sensor  combination,  one  must  maximize  the  utility 
over  the  resource  allocation  into  each  of  the  Nu  sensor  types. 

Let  us  further  examine  the  solution  of  Problem  3  to  see  how  further  simplifi¬ 
cations  can  be  achieved  in  the  resource  allocation  of  the  sensor  combinations  using 
Karush-Kuhn-Tucker  (KKT)  conditions.  Without  loss  of  generality,  we  specifically 
consider  the  case  the  overall  utility  is  the  product  of  the  localization  accuracy  J5]) 
and  lifetime  (JH])  utilities  ( Nu  =  2).  Then,  an  equivalent  problem  is  solve  the  objec¬ 
tive  product  where  A  =  [1  p  0]  ,  and  p  =  a/'y  >  0.  The  overall  utility  is 


K(n)  = 


T 

y  £ji,tht 


T 

y  ui2,tnt 


(17) 


We  start  by  noting  that  the  following  is  a  necessary  KKT  condition  for  the  op¬ 
timal  solution  n*  to  satisfy  over  its  simplex  [Bertsekas  2003;  Epelman  et  al.  2005  : 


hm  >  0, 


dK{h)  dK(n) 

dhm  -  dhi 

nm=n* 


V  l  and  m, 

ni=hy 


(18) 


which  implies  that  at  the  optimum  point,  it  is  not  possible  to  move  resources  from 
any  non-zero  um  to  any  other  zq  without  decreasing  the  utility  K(h).  Using  the 
derivative  of  (ED,  it  is  straightforward  to  see  that  the  necessary  condition  (fTSl)  can 
be  explicitly  written  as  follows: 


T  T 

h*m  >  0,  =>  y  Gmtht  >  y  Guht  V  l  and  m. 

t=i  t-i 


(19) 


where  Gjk  =  6ji:jOJ2,k  +  P&i,k&2,j-  Then,  the  optimal  resource  allocation  solution 
is  given  below  without  a  detailed  derivation: 


1 

0 

(^1  ,j  ~(^l,fc)^’2,fc+p(^2,j  —^2,k)^l,k 

(1  +  P)(^l,i  —&l,k)(&2,j  —&2,k) 


Gjj  >  Gkj’, 

Gkk^_  Gjk 5 

Gjj  G kj  3-11(1  Gkk  ^  Gjk 


(20) 


where  h £  =  1  —  h*.  and  =  0  for  to  ^  j,k.  Hence,  when  Nsc  =  (^)  possible 
combinations  of  the  sensors  in  s  are  considered,  (1201)  can  be  used  to  determine  the 
optimal  resource  allocations.  In  the  general  case,  (fTHl)  can  be  exploited  to  facilitate 
the  resource  allocation  computations. 


3.3.3  Solution  of  Problem  f.  Using  similar  arguments  in  Theorem  1,  the  solu¬ 
tion  to  Problem  4  assigns  resources  to  at  most  Nu  sensor  types.  For  this  problem, 
the  optimization  problem  is  a  linear  program  where  there  are  Nu  linear  constraints 
(Nu  —  1  constraints  are  from  objectives  and  one  cost  constraint).  Therefore,  it  is 
straight  forward  to  solve  Problem  4. 
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Sensor 


Fig.  4.  This  figure  illustrates  the  geometry  of  the  sensor-target  configuration  for  a  monopole 
acoustic  source  moving  along  the  horizontal  direction  with  a  speed  of  v.  The  dashed  lines  represent 
the  acoustic  wave-fronts,  which  create  the  interaction  between  the  target  and  the  sensor. 

Lastly,  once  the  NDS  n*  to  Problems  2,  3,  or  4  is  obtained  via  continuous 
relaxations,  it  is  converted  to  an  integer  solution  via 


ro*  =  floor  ^$— i- 


(21) 


For  large  budgets,  the  NDS  solution  n*  can  be  slightly  different  from  the  solutions 
via  integer  programming  as  the  support  of  the  integer  solution  may  be  greater  than 
Nu.  In  the  integer  solution,  most  of  the  resources  are  used  in  Nu  sensor  types,  but 
a  few  other  sensor  types  help  achieve  a  network  cost  as  close  to  the  budget  $  as 
possible.  In  a  sense,  the  rounding  error  in  (12111  allows  for  a  few  number  of  sensors 
outside  of  the  dominant  Nu  to  appear  in  the  NDS  solution. 

4.  ACOUSTIC  SENSOR  NETWORKS 

In  this  section,  we  focus  on  the  expressions  that  will  lead  to  the  localization  utility 
of  acoustic  sensor  networks  for  position  estimation.  We  first  provide  data  models, 
then  give  expressions  for  the  Fisher  information  matrices  (FIM)  and  the  sensor  re¬ 
ceiver  operating  characteristics  (ROC)  curves  for  position  estimation  using  acoustic 
sensors.  We  assume  that  the  target  and  the  sensors  are  on  the  same  ground  plane, 
measured  in  horizontal  (h)  and  vertical  ( v )  coordinates.  The  target  position  is  de¬ 
noted  as  0  =  \  Oh ,  0V]  and  the  position  of  the  sensor  is  denoted  as  £  =  [  £/,,  ]  . 

The  target  bearing  <f>  is  measured  counterclockwise  with  respect  to  the  sensor’s 
/i-axis  (sensor  orientation:  ip)  and  the  target  range  R  =  \\0  —  £||  is  defined  as  the 
Euclidian  distance  between  the  target  and  the  sensor  (Fig.  [4]). 

Our  results  use  the  analytical  signal  representations  because  time-shifts  in  real 
signals  correspond  to  phase  shifts  in  their  analytical  representations.  We  de¬ 
note  s(f),  x(t),  and  Ui(t)  as  the  complex  envelopes  of  the  source  signal, 

the  source  signal  at  the  array,  the  ith  microphone  additive  noise,  and  the  ith 
microphone  output  signal,  respectively.  To  calculate  target  bearing  or  range,  N 
snapshots  of  the  observed  acoustic  data,  taken  at  times  ti, . . .  ,tjv,  are  used.  We 
note  that  if  the  time  samples  are  sufficiently  apart,  then  successive  samples  of 
the  source  and  the  noise  samples  are  uncorrelated  Davenport  II  and  Root  1958| 
IWeiss  and  Weinstein  1983],  We  model  the  source  signal  samples  as  i.i.d.,  zero  mean, 
complex  circularly  symmetric  Gaussian  random  variables  with  variance  a2  and  the 
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noise  samples  with  i.i.d.,  zero  mean,  complex  circularly  symmetric  Gaussian  ran¬ 
dom  variables  with  covariance  a2 1. 


4.1  Range  Sensors 

We  first  discuss  the  range  estimation  of  a  narrow-band  source  using  an  omni¬ 
directional  microphone  in  an  isotropic  medium.  We  assume  that  there  are  no 
multipath  effects.  Assuming  spherical  propagation,  we  write  the  envelope  of  the 
microphone  output  signal  at  the  target  narrow  band  frequency  fw  as  follows, 
|Morse  and  Ingard  1968[  |Johnson  and  Dudgeon  1993  : 

y(t)  =  x(t )  +  n{t)  =  e~J  +  n(t),  (22) 

where  j3(t)  =  1  +  ^  cos  ip  is  called  the  Doppler  shift  factor  and  c  is  the  speed  of 
sound.  Based  on  our  signal  and  the  noise  signal  assumptions,  it  is  straightforward  to 
prove  that  the  microphone  output  signal  y{t)  also  has  an  i.i.d.  zero  mean  circularly 
complex  Gaussian  distribution  with  variance  a2  =  a2  +  a2,  where  er^  =  Now, 
we  denote  the  JV-sample  root-mean-squared  microphone  output  as  e: 


\ 


N 


i>(*i)  i2= 


E 


(  Prealfc 


r?  V  ay/2 


+ 


y ?  (ti) 


5/2 


V2N 


z, 


(23) 


where  we  define  z  as  the  second  square-root  summation  term  in  (1231) .  The  variable 
z  has  a  Chi  distribution  pz(z)  with  2N  degrees  of  freedom  (Evans  et  al.  200f)|. 

Although  we  now  have  an  analytical  expression  for  e,  we  apply  the  Laplacian  ap¬ 
proximation  to  further  facilitate  the  derivations  of  the  ROC  curves.  The  Laplacian 
approximation  uses  the  mode  and  the  Hessian  of  the  log  likelihood  at  the  mode  to 
approximate  it  with  a  Gaussian.  The  resulting  approximation  is  (IV  ^>1) 


£ 


p(£)  =  U 


■  N 


+  a2 


,  2 
a.x  +  a 


AN 


(24) 


The  Fisher  information  matrix  F  is  defined  as  follows  [Lehmann  and  Casella  1998]: 


F(6)  =  j  p{e\6) 


d\ogp{e\0) 


80 


<91ogp(e|0) 


80 


de. 


(25) 


We  use  p(e)  in  (l24l)  to  determine  an  approximate  FIM  of  the  target  position  estimate 
given  below  without  a  detailed  derivation: 


Fi(0)  = 


AN(a2/a2f 
R2  (a2 /a2  +  R2)2 


cos  t[> 
sin  tjj 


x  [  cos  tp  sin  ip  1  , 


(26) 


where  %p  is  the  sensor  bearing  with  respect  to  the  target  (see  Fig.  [4])  and  it  is 
assumed  that  a//3  w  1. 

In  the  case  of  wideband  sources,  if  the  observation  period  of  the  N  samples  is 
much  larger  than  the  inverse  bandwidth  of  the  source  signal  s(f),  but  short  enough 
so  that  the  target  is  approximately  stationary,  then  the  discrete  Fourier  transform 
(DFT)  coefficients  of  the  signal  are  statistically  uncorrelated,  [Davenport  II  and  Root  1958 
IWeiss  and  Weinstein  1983;  IChow  and  Schulthciss  1981].  Using  the  properties  of 
the  FIM,  [Van  Trees  1968|.  the  wideband  FIM  is  a  summation  of  the  individual 
narrow-band  FIM’s. 
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To  derive  the  receiver  operating  characteristics  (ROC)  curves  of  the  range  sensors, 
we  start  with  (PHI).  Note  that  when  N  3>  1,  we  can  write  e  «  ay  +  -^=A/r( 0, 1)  ~ 

A/~(0,1) 

aye  yuv  .  Hence,  a  multiplicative  noise  model  on  the  envelope  is  more  appropriate 
than  an  additive  model.  Denoting  £  =  loge,  we  have 

£~Af(\og<7v,-±^  (27) 

with  the  following  hypotheses: 


Hq'-  Gy  —  Gni 

Hi:  Gy  =  \j  G2  +  G2  >  Gn. 


(28) 


In  this  case,  likelihood  ratio  test  is  equivalent  to  a  simple  linear  threshold  detector: 


Hi 

£  ^  V,  (29) 

Ho 

and  this  detector  is  also  uniformly  most  powerful  (UMP)  [Van  Trees  19681  iLehmann  and  Casella  1998  . 
It  is  easy  to  verify  that  the  ROC  curve  is  determined  by 

Pd  =  $  (V1  (Pf)  -VN(\og(G2  +  g2x)  -  log  a2))  ,  (30) 

where  $(•)  is  one  minus  the  cumulative  distribution  function  of  the  N  (0,  l)-random 
variable,  Pd  and  Pf  are  the  detection  and  false  alarm  probabilities,  respectively. 


4.2  Bearing  Sensors 

Bearing  sensors  determine  the  direction-of-arrival  (DOA)  of  a  target  by  exploiting 
the  time-delay  information  present  in  the  elements  of  the  acoustic  array  [Johnson  and  Dudgeon  1993 
When  the  number  of  acoustic  data  samples  used  to  calculate  the  bearing  are 
sufficiently  high,  the  estimated  bearing  has  the  following  Gaussian  distribution, 

[Bell  et  al.  19961  ILiu  et  al71|  Johnson  and  Dudgeon  1993|: 

<t>  =  tan-1  ^  _  ^)  +  n<A:  (31) 

where  the  bearing  noise  n $  depends  on  the  actual  microphone  noise  on  the  array 
elements. 

To  determine  the  relationship  between  the  bearing  noise  to  the  additive  mi¬ 
crophone  noise,  consider  a  narrow-band  plane  wave  signal  impinging  on  an  m- 
microphone  planar  acoustic  array.  We  first  derive  the  Fisher  information  ma¬ 
trix  for  the  bearing  cj>  and  approximate  the  variance  of  using  the  inverse  of 
theFIM,  |  Johnson  and  Dudgeon  199~3|IStoica  and  Nehorai  19891  IStoica  et  al.  200H 
iBell  et  al.  1996],  The  wavenumber  vector  k  is  given  by 

k  =  ^  w  u,  u  =  T  cos  4>  sin (/>]',  (32) 

c  L  J 

where  fw  is  the  narrow-band  source  temporal  frequency,  c  is  the  speed  of  propa¬ 
gation,  and  u  is  the  Cartesian  bearing.  The  microphone  locations  in  the  acoustic 
array  are  given  by  the  matrix  D  defined  as  follows: 


D  —  (  dl  .  .  .  drn  ] 


dh,l  •  •  •  dh:n 


dii 


d/i 


(33) 
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where  di  has  the  location  of  the  ith  microphone  (i  =  1, ...  ,  to)  in  the  Cartesian 
coordinate  system.  The  time  delay  at  the  itli  microphone  relative  to  the  origin  of 
the  array  is  given  by  Tj  =  u  .  Then,  the  array  output  can  be  written  as 

y(t)  =  a(u)x(t)  +  n(t),  (34) 


where  y{t)  =  [  j/i  (t)  ...  ym(t)  ]',  n(t)  =  [  m(t)  ...  nm(t)  ] a(t)  =  [  e-j27r/”T1  ...  e 
Under  stationarity  and  Gaussian  assumptions,  it  can  be  shown  that  the  Fisher  in¬ 
formation  matrix  for  a  narrow-band  source  is  given  by  [Bell  et  al.  1996] : 
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where  the  matrix  term  involving  the  position  matrix  can  be  shown  equal  to  m(ar /2)I 
for  uniform  circular  arrays  where  a  is  the  array  radius.  Then,  it  is  straightforward 
to  determine  the  Fisher  information  matrix  for  the  bearing  </>  via  the  following 
transformation,  e.g.,  see  [Van  Trees  1968]: 

F^)  =  (V^uf  F(u)  (V^w) .  (36) 


Note  that  when  the  acoustic  array  has  a  uniform  circular  geometry,  does 

not  depend  on  <f>.  Using  m  and  (1361) .  we  present  the  array  FIM’s  for  the  target 
position  without  a  detailed  derivation  (see  Fig.  [4]): 


Fm(  0) 
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sin  ip 
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(37) 


The  wideband  version  for  the  bearing  sensors  is  a  summation  of  the  corresponding 
narrowband  FIM’s  similar  to  the  range  sensors. 

To  derive  the  ROC  curves  for  bearing  sensors,  we  use  the  following  detector: 


Hi 


maxfra 

m 


^  v , 

Ho 


(38) 


where  £m  =  log£m  for  the  mth  sensor.  Probability  that  the  maximum  of  m  statisti¬ 
cally  independent  Gaussian  random  variables  £m  with  the  mean  log  er2  and  variance 
^  exceeds  the  threshold  rj  is  given  by 

m 

Pf  =  1  -  [Q  [l  (V4N(y-  0.5 logo-2))]  .  (39) 

k= 1  L 

Similarly,  the  detection  probability  is  given  by 

Pd  =  1  -  [l  -  Q  (Q-1  ((1  -Pf)jr)-y/N  (log  (a2  +  a2x)  -  log  a2))]  ™  .  (40) 


5.  SIMULATIONS 

5.1  Example  Dynamic  Programming  and  Approximate  Branch-and-Bound  Solutions 

Using  a  synthetic  example,  we  demonstrate  the  dynamic  programming  solution 
of  the  NDS  problem  and  also  compare  them  with  the  approximate  solutions  to 
establish  their  closeness.  This  comparison  is  quite  important  as  the  approximate 
solutions  obtained  by  continuous  relaxations  sometimes  fail  to  be  useful  in  finding 
a  branch-and-bound  solution  |Nemhauser  and  Wolsey  1988| .  We  use  the  following 
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synthetic  problem:  T  =  6,  ct  =  1  +  t,  ft  =  t2,  S  =  1.2,  7  =  2.2,  $  =  500,  and 
R*  =  [  1,  2,  2,  2,  3,  3  ]'.  With  these  parameters,  sensor  type  6  has  the  maximum 
whereas  sensor  type  5  has  the  maximum 

Figures  Ha)  and  (b)  show  the  evolution  of  the  dynamic  programming  solution 
for  Problem  1  and  the  Pareto  frontier  surface.  In  Fig.  03a),  u±  is  maximized  at 
N  =  72  where  n  =  [0  1  0  0  0  71]'  (Problem  2).  When  N  <  71,  the  optimal 
n  consists  of  only  s6  (see  Fig.  [5](c)  (bottom)).  Due  to  the  integer  nature  of  the 
problem,  N  =  72  gives  slightly  better  result,  because,  when  N  =  71,  the  cost 
c  =  71  x  7  =  497  <  $  results  in  l  =  71  x  62  =  2556,  whereas,  when  TV  =  71  + 1  =  72, 
c  =  71x7  +  lx3  =  500  <  $  results  in  l  =  71  x  62  +  1  x  22  =  2560.  Note 
that  this  is  less  than  l  =  (71  +  3/7)  x  62  =  2571.4  if  it  were  possible  to  buy 
fractional  sensors.  However,  most  resources  are  only  spent  on  the  best  performing 
sensors  for  their  cost  as  expected.  In  Fig.  Eta),  M3  is  maximized  at  N  =  84  where 
n  =  [  1  0  0  0  83  0]'.  As  expected,  most  of  the  resources  are  spent  on  sensor  type 
5.  Figure  Hc)  demonstrates  the  distribution  of  the  resources  (top)  at  a  constant 
coverage  V  while  trading  lifetime  L  for  localization  utility  U  starting  from  maximum 
utility,  (middle)  at  a  constant  lifetime  while  trading  coverage  V  for  utility  U  starting 
from  the  maximum  coverage,  and  (bottom)  at  the  boundary  of  the  T  solutions  where 
the  lifetime  L  is  traded  for  coverage  V.  In  all  cases,  sensor  type  3  rarely  gets  assigned 
any  resources.  This  is  because  the  rows  of  Cl  are  locally  convex  around  sensor  type 

3,  hence  sensor  type  3  gets  dominated  by  the  combinations  of  sensor  types  2  and 

4.  In  fact,  on  the  whole  Pareto  frontier,  n%  <  1,  where  it  is  only  assigned  resources 
due  to  the  integer  nature  of  the  problem.  Hence,  for  all  practical  purposes,  sensor 
type  3  could  have  been  eliminated  before  any  solution  is  attempted. 

For  the  solution  of  Problem  3,  it  is  straightforward  to  see  that  the  convexity 
conditions  for  the  approximate  solutions  are  satisfied.  Analytical  formula  (l20l)  for 
the  solution  of  Problem  3  results  in  iq  =  0.1321  and  u§  =  0.8679  with  um  =  0  for 
m  =  2, ...  ,5.  The  continuous  resource  variables  um’ s  correspond  to  n  estimates  as 
n\  =  33.0317,  uq  =  61.9910,  and  nm  =  0  for  m  =  2, . . . ,  5  via  nm  =  um%/cm.  The 
dynamic  programming  solution  in  Algorithm [T]  results  in  n*  =  [  33  0  0  0  0  62]', 
which  closely  exhibits  the  boundary  effect  phenomena  as  predicted  by  the  branch- 
and-bound  solution. 

Figures  m  a)  and  (b)  illustrate  the  evolution  of  the  dynamic  programming  so¬ 
lution  and  Fig.  He)  plots  the  Pareto  frontier  for  the  solution  of  Problem  3.  In 
Fig.  Hb),  the  solution  of  Problem  3  is  shown,  where  N  =  95,  where  we  decrease 
the  maximum  utility  by  1.305  times  to  increase  the  lifetime  by  1.418.  In  Fig.  He), 
we  show  other  operating  points  on  the  Pareto  frontier.  For  example,  to  increase  the 
lifetime  3.47  times,  the  utility  needs  to  be  decreased  10  times.  Finally,  we  observed 
that  all  the  points  on  the  Pareto  frontier  exhibit  the  boundary  behavior  similar  to 
what  we  expected  for  the  solution  of  Problem  3. 

5.2  Acoustic  Sensor  Network  Design 

In  this  section,  we  use  our  results  to  design  an  acoustic  sensor  network  for  position 
estimation.  Typical  target  acoustic  time-frequency  signatures  are  shown  in  Fig.  0 
Based  on  the  typical  vehicle  time-frequency  characteristics,  we  model  the  power  of 
a  vehicle  as  wideband  with  flat  spectral  characteristics  between  0  —  300Hz.  The 
intersensor  spacing  of  the  uniform  circular  acoustic  arrays  is  determined  using  the 
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N  @  V  on  Upper  Boundary 


(c) 


Fig.  5.  (a)  The  dynamic  programming  solution  via  Algorithm  [3  for  Nu  =  3.  (b)  Pareto  efficient 
frontier  is  shown  with  different  operating  points.  The  utility  U  and  the  lifetime  L  are  plotted  in 
logarithmic  scale  and  are  normalized.  The  maximum  U  is  —  1.  Also,  L(Nfiin in)  is  1.  (c)  Resource 
distribution  is  shown  for  various  paths  on  (a).  Darker  colors  imply  more  resources.  Note  that 
sensor  type  3  is  rarely  assigned  resources. 


Fig.  6.  (a)  The  dynamic  programming  solution  via  Algorithm  [T]  for  Nu  =  2.  (b)  Corresponding 
solution  of  Problem  3  is  illustrated  (c)  Pareto  efficient  frontier  is  shown  with  different  operating 
points.  The  dashed  line  shows  the  supporting  line  with  the  maximum  slope  going  through  the 
origin. 


bandwidth  of  the  source  signals.  We  constrain  the  sensor  sizes  to  be  less  than  or 
equal  to  lm  for  practical  reasons. 
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(a)  Compact  (b)  SUV 


Time  in  [s] 

(c)  Truck 


Fig.  7.  Time  frequency  images  in  dB  scale  for  (a)  a  compact  vehicle  (b)  a  sports  utility  vehicle 
(SUV),  and  (c)  a  truck,  each  moving  at  approximately  40km/h.  The  sampling  rate  for  the  micro¬ 
phone  is  Fs  =  44100Hz.  The  microphone  has  flat  spectral  response  up  to  1000Hz.  A  notch  filter 
is  used  to  minimize  the  electrical  interference  at  60Hz.  At  moderate  vehicle  speeds,  most  of  the 
vehicle  signal  energy  lies  within  0-300Hz  region.  Moreover,  most  of  the  energy  is  approximately 
uniformly  distributed  in  addition  to  the  concentration  around  engine  cylinder  firing  frequencies. 


To  characterize  the  typical  vehicle  signal  power  levels,  we  use  the  published  mea¬ 
surements  of  the  U.S.  Department  of  Transportation  Federal  Highway  Administra¬ 
tion  |Department-of- Transportation  |.  For  vehicles,  a  functional  form  of  the  vehicle 
noise  power  is  given  by 

of  =  10  log10  (lO(C/10)  WA/10)10(S/10))  +  10  log10  (152)  ,  (41) 

where  v  is  the  vehicle  speed  in  units  of  km-per-hour,  A  =  41.74,  B  =  1.149,  and 
C  =  50.13.  In  (ITTjl.  the  second  logarithmic  term  accounts  for  the  measurement 
distance,  which  is  15m.  Hence,  m  illustrates  the  average  vehicle  signal  power 
for  automobiles,  medium  trucks,  heavy  trucks,  buses,  and  motorcycles.  The  signal 
power  has  a  flat  level  at  <r 2  ss  73dB  when  v  <  lOkrn/h  due  to  the  engine  and  ex¬ 
haust,  then  has  a  linear  increase  when  v  >  lOkm/h  (of  ss  107dB  at  v  =  lOOkm/h) 
attributed  to  the  tire  and  pavement  noise.  For  our  example  design  in  this  section, 
we  exclude  the  heavy  trucks,  e.g.,  18- wheelers,  and  assume  that  the  average  sig¬ 
nal  power  over  the  ambient  noise  is  uniformly  distributed  between  50-80dB  (i.e. , 
101og10  (of/cr2)  ~  U  (50,80)). 

We  use  ci  =  5  and  ct  =  25  +  5 1  to  design  an  acoustic  sensor  network  using  a 
budget  of  10K  units  to  cover  a  surveillance  region  of  A  =  1500mx  1500m  and  impose 
that  the  connectivity  probability  pc  of  the  network  to  be  greater  than  0.99.  In 
cost  schedule,  acoustic  arrays  are  more  expensive  than  range  only  sensors  because 
(i)  they  require  special  deployment  mechanisms  and  (ii)  they  require  additional 
orientation  calibration.  It  is  easy  to  see  that  the  convexity  assumption  is  satisfied 
for  u>2,t- 

To  determine  Qij,  we  first  estimated  a  target’s  position  situated  at  the  center  of 
a  disk  of  radius  maxR^  over  1000  Monte  Carlo  realizations  to  determine  estimation 
variances  of  and  er2  for  each  sensor  type  by  varying  nt .  To  estimate  the  position 
variances  for  range  only  sensors,  we  generated  acoustic  data  using  (12211  with  v  =  0, 
used  range  estimates  from  the  sensors  with  positive  detections  according  to  (l29l). 
then  used  a  Newton-Raphson  recursion  with  the  known  sensor  positions  and  esti¬ 
mated  a  target  position.  We  used  R\  =  50m  as  the  operational  range  of  the  range 
only  sensors.  To  estimate  the  position  variances  for  bearing  only  sensors,  we  used 
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the  inverse  of  the  wideband  version  of  the  FUV@  in  (1551)  to  generate  bearing  esti¬ 
mates  for  each  array,  a  Newton-Raphson  recursion  with  the  known  sensor  positions 
and  estimated  a  target  position  (R%  =  300m).  We  then  took  the  inverse  of  af  -f  <r^ 
for  each  sensor  type  to  estimate  the  ft  and  7  via  (11- 

Figures  Ufa)- (f)  show  some  system  simulations  using  the  single  sensor  and  pair¬ 
wise  sensor  combinations  (solid  curves).  In  the  figures,  the  dashed  lines  represent 
the  results  of  our  approximation  in  © ,  where  we  have  estimated  7  =  2.18420  The 
corresponding  ft  estimates  are  shown  in  Fig.  ©a)  and  (b).  In  Fig.  ©c),  we  plot 
the  performance  per  cost  curve  u>i =  ft/ct  for  the  assumed  cost  structure.  Unfor¬ 
tunately,  wM  does  not  have  a  convex  form,  however,  it  can  be  closely  upper  and 
lower  bounded  by  straight  lines  as  marked  with  stars  and  squares  in  Fig.  ©c). 


Fig.  8.  System  simulations  (solid  lines)  and  the  solution  of  (dashed  lines),  (a)  s  1  alone,  (b) 
S5  alone,  (c)  sio  alone,  (d)  si  and  S2 ■  (e)  S3  and  sq.  (f)  s  1  and  s n. 


Figures  [TOT a)  and  (b)  illustrate  the  evolution  of  the  dynamic  programming  solu¬ 
tion  and  Fig.  mac)  plots  the  Pareto  frontier  for  the  acoustic  sensor  network  design 
along  with  the  solution  of  Problem  3  and  the  corresponding  upper  and  lower 
bounds  using  (l20l).  Figure Q©d)  shows  the  resource  allocation  on  the  Pareto  fron¬ 
tier.  The  dynamic  programming  solution  of  Problem  3  results  in  m  =  720, 
rin  =  80  and  zero  for  the  rest  of  the  sensor  types.  These  numbers  correspond 


3  The  wideband  version  uses  a  flat  spectrum  assumption  and  is  a  summation  of  (I36|)  at  each  FFT 
frequency  determined  by  the  sampling  rate.  This  FIM  is  known  to  be  tight  for  bearing  errors 
when  the  number  of  data  samples  is  high,  which  is  the  case  here.  For  example,  MUSIC  algorithm 
achieves  this  bound,  e.g.,  see  [Stoica  and  Nehorai  1989]. 

4The  value  of  the  parameter  7  depends  on  the  target  SNR  range.  When  the  SNR  decreases,  7 
increases  and  vice  versa. 
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to  approximately  2.5-Si  and  10.05-Sn  sensors  in  areas  of  size  t:R\2  and  irR^2, 
respectively.  Hence,  each  point  on  A  is  well  covered  at  least  by  3  sensors.  The 
minimum  radio  transmission  range  should  be  chosen  as  rj*ran  =  100.51m  via  0 
to  guarantee  a  connectivity  probability  of  pc  =  0.99  (d  =  1).  The  lower  bound 
results  in  u\  =  0.3549  and  un  =  0.6451  via  (fl9l)  corresponding  to  n\  =  709.7365 
and  Tin  =  80.6415.  If  the  resources  for  0.6415-Sn  sensors  are  spent  on  sensors  of 
si  type,  the  approximate  design  results  in  ni  =  720  and  nn  =  80,  which  inciden¬ 
tally  coincides  with  the  dynamic  programming  solution.  The  upper  bound  results 
in  m  =  0.4098  and  «n  =  0.5902  via  (fl9ll  corresponding  to  ni  =  819.5500  and 
nn  =  73.7781. 


(a)  (b)  (c) 


Fig.  9.  (a)  It  is  instructive  to  plot  /T  as  it  demonstrates  the  relative  performance  of  the  sensor 
networks  that  consist  of  only  one  sensor  type  sm.  There  is  approximately  a  cubic  increase  up  to 
sn,  where  the  arrays  hit  the  radius  limit  of  lm.  After  sn,  the  increase  is  approximately  m1A. 
The  cubic  increase  can  be  explained  by  (i)  the  aperture  gain  due  to  the  a2-term  in  the  FIMs,  and 
(ii)  SNR  gain  in  the  FIMs.  At  any  single  SNR,  only  the  aperture  gain  is  present,  (b)  Inherent 
performance  /m  is  shown,  (c)  Performance  per  cost  u)ift  for  the  given  cost  structure  is  shown. 
Convex  upper  and  lower  bounds  are  marked  with  squares  and  stars.  The  dynamic  programming 
solution  uses  the  10  times  the  values  of  the  dots,  which  are  integer,  that  approximate  Acm. 


6.  CONCLUSIONS 

In  summary,  we  presented  mathematical  programming  and  approximation  algo¬ 
rithms  for  the  NDS  problem  for  acoustic  sensor  networks  using  position  estimation 
as  an  example.  We  provided  approximation  algorithms  that  can  bound  the  results 
as  well  as  methods  to  identify  dominated  sensors  to  alleviate  computation.  For  a 
given  deployment  region  size,  our  dynamic  programming  solution  can  calculate  the 
exact  integer  Pareto  frontier  of  the  sensor  network  utility  at  the  desired  probabil¬ 
ities  for  d-connectivity  and  fc-coverage.  However,  the  continuous  relaxations  are 
shown  to  be  quite  useful  in  obtaining  quick  solutions  to  the  problem. 

To  the  Pareto  problems  in  this  paper,  other  linear  matrix  inequalities  can  be 
added  without  changing  the  solution  strategies.  As  an  example,  the  total  weight  of 
the  sensor  network  might  also  be  constrained  in  space  exploration,  where  a  mobile 
robot  carries  the  sensors  for  deployment.  This  results  in  an  affine  constraint,  which 
can  be  treated  as  convex  or  concave  as  required. 

As  future  work,  we  will  extend  the  formulation  in  this  paper  to  hypothesis  test¬ 
ing.  Suppose  that  based  on  the  state-of-nature,  the  binary  classification  problem  is 
linear,  where  we  choose  Hq,  when  bT 0  >  0,  and  Hi  otherwise  |Duda  et  al.  2001], 
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Fig.  10.  (a)  and  (b)  Dynamic  programming  solution  of  the  NDS  problem  with  the  minimum  cost 
are  plotted,  (c)  Pareto  efficient  frontier  is  shown  with  different  operating  points.  Due  to  the 
non-convexity  of  not  all  the  point  on  the  Pareto  frontier  can  be  determined  by  scalarization. 
The  obtained  upper  and  lower  bounds  (the  square  and  the  star,  respectively)  are  quite  close  to 
the  dynamic  programming  solution,  (d)  The  resource  distribution  on  the  Pareto  frontier  is  shown. 


In  this  case,  it  can  be  shown  that  the  Bayesian  design  criteria  becomes  the  fol¬ 
lowing  [Toman  1996] :  U(n)  oc  — tr  |bbTASr~1(n)|.  For  sensor  networks  with 
acoustic  and  video  modalities,  the  discriminant  features  can  be  the  acoustic  am¬ 
plitude  and  the  object  size  to  classify  vehicles,  e.g.,  into  compact  vs.  SUV  cate¬ 
gories  [Cevher  et,  al.  20071  ICevher  et,  al.  2007] .  We  plan  to  investigate  Pareto  fron¬ 
tiers  for  joint  parameter  estimation  and  classification.  The  resulting  NDS  problems 
will  then  address  to  a  broader  class  of  sensor  network  problems  and  present  further 
intellectual  challenges  for  mathematical  programming. 
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